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Abstract 

We address the study of a class of ID nonlocal conservation laws from a numerical point 
of view. First, we present an algorithm to numerically integrate them and prove its 
convergence. Then, we use this algorithm to investigate various analytical properties, 
obtaining evidence that usual properties of standard conservation laws fail in the nonlocal 
setting. Moreover, on the basis of our numerical integrations, we are lead to conjecture 
the convergence of the nonlocal equation to the local ones, although no analytical results 
are, to our knowledge, available in this context. 
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1 Introduction 

Conservation laws with nonlocal fluxes have appeared recently in the literature, arising nat- 
urally in many fields of application, such as in crowd dynamics (see O [H O [15] and the 
references therein) , or in models inspired from biology, see [El El El EOl • 

In this paper, we initiate the study of these equations from a numerical point of view. 
First, we prove the convergence of a finite volume algorithm to numerically integrate a class 
of one-dimensional conservation laws with a nonlocal flow. Then, we use this algorithm to 
show peculiar properties of these nonlocal equations and, in particular, how they differ from 
the usual local ones. 

Consider the scalar equation 

I ^fn ^ ^^o?' f ' * ''^^ ^ ° (t, X) G M+ X R (1.1) 

1 p(0, X) = p^[x) 

which slightly extends, in the ID case, the class of equations considered in [4^^. We present 
below a numerical scheme to integrate and prove its convergence. As a byproduct, we 

also establish an existence result for fjl.ll) , thus slightly extending [H Theorem 2.2] in the ID 
case. 

This numerical algorithm is then implemented and used to investigate various properties 
of (II. ip . First, we provide evidence that the usual Maximum Principle for scalar conservation 
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laws fails in the case of (II. ip . Another integration shows that the total variation of the solution 
to (jl.ip may well sharply increase, contrary to what happens in the standard local situation. 
Remark that both these examples are in agreement with the estimates we rigorously obtain 
on the approximate solutions. 

Of particular interest is the limit 77 ^ 5, 5 being the Dirac measure centered at the origin. 
Numerical integrations show that the solutions to (II. ip converge to that of 

^' "^^^^ ^ ° (t, X) G M+ X R (1.2) 
p(0, x) = p^\x) 

although no rigorous proof of this convergence is, to our knowledge, known. Remark that in 
the nonlocal case, well posedness results are available also in the case of systems in several 
space dimensions, see [5l[6]. Hence, the ability of passing to the limit r] ^ 6 might help in the 
search for analytical results about systems of conservation laws in several space dimensions. 

Let us make the following remarks. First, the scheme presented below has an associated 
CFL condition. The CFL condition is often interpreted through a comparison between the 
numerical propagation speed and the analytical one, see for instance [131 § 4.4, p. 68]. In the 
present nonlocal case (II. ip , information propagates at an infinite speed, due to the presence 
of the term 77 * p. Nevertheless, also in the nonlocal case (|l.ll) a suitable CFL condition plays 
a key role, see (12. 4p . 

Second, the scheme presented below is not monotone in the sense of the usual defini- 
tion jT3l Formula (12.42)], as follows from the integration in Section 1X21 There, both constant 
initial data p = and p — 1 yield constant solutions, but the initial datum (j33j), although 
it attains values in [0, 1], yields a solution exceeding 1. Nevertheless, the scheme (12. 6p enjoys 
several properties of monotone schemes, proved in the lemmas in Section [2l 

Remark 1.1. Throughout this work, we follow the usual habit of referring to (II. ip as to 
a nonlocal equation and, hence, to the standard case (II. 2p as to the local case. However, 
whenever the support of 77 is bounded, it might seem more appropriate to call (II. ip a local 
equation and (II. 2p the punctual case. 

The next section deals with the definition of the algorithm and with the statement of 
the estimates which ensure its convergence, as well as the entropicity of the limit solution. 
Section E] deals with various numerical integrations of (jl.ip . All proofs are deferred to the 
last Section m 



2 Main Results 

Throughout, we set = [0, +oc[. 

As a starting point, we state what we mean by solution to (II. ip , see also [4^ Definition 2.1]. 

Definition 2.1. Let T > 0. Fix G L^(R;R). A weak entropy solution to ([Ll]) on [0,r] 
is a hounded measurable Kruzkov solution p G ([0, T]; L|^^(R; R)) to 
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For the definition of Kruzkov solution, see for instance [71 Paragraph 6.2] or \T2, Definition 1]. 
Here, as usual, 

* ry) (x) = / . 

Jr 

Remark that the assumptions 



fee' 



X K X 



^; e (C^nw 



l.oo 



and 



; R) and 



sup < +OC 

t,x,p 

sup < C\p\ 

t,x 

dlj{t,x,p) < c\p\ 



sup 

t.x 



\/t,x f{t,x,0) 
rye (C2nW2'°° 



(2.1) 



(2.2) 



ensure that the transport equation in Definition 12.11 fits in Kruzkov framework, see [71 [T^ . 
From the modeling point of view, it is natural to require that the kernel 77 attains only positive 
(or non-negative) values. However, this requirement is not necessary for the analytical results 
below. 

Below, Remark 12.31 and Lemma 12.51 provide uniform bounds on the solution to (|l.ip 
under conditions (I2.ip - fl2.2p on the equations and for data in L^. Therefore, the apparently 
strong requirement ||^p/||loo < +^ can be easily relaxed to 



sup 

teR+ ,xeR, pe[-M,M] 



\dpf{t,x,p)\ < +00 



for a suitable positive M. Moreover, the usual sublinearity condition sup^^, \dxf{t,x,p)\ < 
+ \p\) takes the form sup^ 2, \dxf(t^ p)\ < C\p\ in (12. ip due to the assumption /(t, x, 0) = 
for all t and x. 

Introduce a uniform mesh with size h along the x axis and size r along the t axis. Through- 
out, we assume that 

h < 1/C 

with C as in (12. ip and that the following CFL condition is satisfied: 



(2.3) 



^ (l + 2||9p/|| 



Loo 



< 



1 



(2.4) 



where, as usual, X — r/h. 

Consider the following Lax-Friedrichs type scheme: 



P] = T / P"(x)dx 

where the numerical flux fjYi/2 (ES]) is given by 



(2.5) 



-^(P2-Pl) 



(2.6) 
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Here, the convolution is computed through a standard quadrature formula using the same 
space mesh, as follows 

= 5^ ^ Pk+1/2 Vj+l/2-k (2.7) 

where ^ suitable convex combination of p^, p^^^ and Vj+1/2 — j[ Jx^^^ r]{x)dx, for 

instance. 

The next three lemmas provide the basic properties of the algorithm (12. 5p , namely posi- 
tivity, and bounds. All proofs are deferred to Section [H 

Lemma 2.2 (Positivity) . Let conditions (12. II) - (12. 2 p hold. Assume that h and r satisfy (12. 3p 
and the CFL condition (12. 4p . If > for all then the approximate solution constructed 
by the algorithm (12. 5p is such that p^ > for all j and n. 

Remark 2.3. The proof of the above lemma clearly shows that if we assume p^ < 0, then < 
for all n. Moreover, under the same assumptions (I2.ip - (l2.2p - (l2.3p - (l2.4p , a straightforward 
modification of the proof of Lemma 12.21 ensures that if there exists a p G such that 
/(t, X, p) — 0, then the inequality p^ > p, respectively p^ < p, implies that > p, respectively 
P^ < p, for all n. 

Lemma 2.4 (L-*^ bound). Let conditions (I2.ip - (l2.2p hold. Assume that h and r satisfy (12. 3p 
and the CFL condition (12.41) . If Pj > for all j, then the approximate solution constructed 
by the algorithm (12. 5p satisfies 

IIp"IIli<IIp°IIli- 



Lemma 2.5 (L^ bound). Let conditions (I2.ip - (l2.2p hold. Assume that h and r satisfy (12. 3p 
and the CFL condition (12. 4p . If Pj > for all j, then the solution constructed by the algo- 
rithm (12. 5p satisfies 

\\P IIloo ^ \\p IIloo e , 

where C depends on C in 1^2. on various norms of f^v^rj and on the norm of the initial 
datum, see (14. 5p . 



The next result concerns the bound on the total variation of the approximate solution 
constructed in (12. 5p . In the standard Kruzkov case, when the flow is independent from t and 
X, the total variation of the solution is well know to be a non-increasing function of time, 
see [U Theorem 6.1]. Here, on the contrary, the total variation and the L^ norm of the 
solution to (jl.ll) may well sharply increase due to the nonlocal terms, even when the flow is 
independent from t and x, see Section | 



Proposition 2.6 (Total variation bound). Let conditions (j2.1l) - (j2.2l) hold. Assume that h 
and T satisfy (|2.3I) and the CFL condition (|2.4I) . If pj > for all j, then the approximate 
solution constructed by the algorithm (|2.5I) satisfies the following total variation estimate, for 
all n > 0; 




e^'^ , (2.^ 



where the constants /Ci and /C2 depend on C in 112. on various norms of f,v,r] and of the 
initial datum, see (|4.18l) . 



4 



A first consequence of the bound on the total variation is the L^-Lipschitz continuity in 
time of the approximate solution, proved in the following lemma. 

Lemma 2.7 (L-'^-Lipschitz continuity in time). Fix a positive T. Let conditions (|2.1l) - (|2.2p 
hold. Assume that h and r satisfy (12. 3p and the CFL condition (|2.4I) . If pj > for all then 
the approximate solution constructed by the algorithm (12. 5p is an h^-Lipschitz continuous 
function of time, in the sense that for any n, m G N such that nr <T and mr <T , 

IIp^-P^IIli <C(T)\n-m\T 

where the quantity C{T) grows exponentially in time and depends on C in 112. on various 
norms of f^v^rj and of the initial datum, see (14.191) . 



The bound proved in Lemma [2. 5^ the total variation bound proved in Proposition 12.61 
and the uniform continuity in time that follows from Lemma lZTl allow to apply Helly Theorem, 
for instance in the form of [ll Theorem 2.6], to the sequence of approximate solutions con- 
structed through (|2.5I) . A straightforward limiting procedure, see for instance [1, Section 6.2], 
thus ensures the existence of weak solutions to the Cauchy problem for (jl.ll) . 

To obtain uniqueness, we prove that the approximate solutions (j2.5l) also satisfy a discrete 
entropy condition. To this end, define for each /c G M the Kruzkov numerical entropy flux as 

Fj+l/2ipl^P2) = fjVl/2(Pl V fc, P2 V fc) - fjVi/2(pi Ak,p2Ak), (2.9) 
where a\/ h — max(a, h) and a /\h — min(a, h). 

Proposition 2.8 (Discrete entropy condition). Let conditions (|2.ip - (|2.2l) hold. Assume that 
h and r satisfy (j2.3l) and the CFL condition (j2.4l) . If > for all j, then the approximate 



solution constructed by the algorithm (|2.5I) verifies the discrete entropy inequality 



(2.10) 



for all k G 



3 Numerical Integrations 
3.1 A Nonlocal Traffic Model 

The classical Lighthill-Whitham [14j and Richards (LWR) model for vehicular traffic 
consists of the continuity equation dtp + dx{pV) = supplied with a suitable speed law 
V = V(p). Here, as usual, t is time, x an abscissa along a rectilinear road with neither entries 
nor exits and p G [0, 1] is the (average) vehicular density. 
Equation (jl.ll) with 

IT /O 

f{p) = P(l - P) ^(r) = y^ax (1 - r) and r]{x) = a {{x - a){b - x)) X[^^^j(^) , (3.1) 

where l^max > 0, can be used as an LWR-type macroscopic model for vehicular traffic, where 
drivers adjust their speed according to the local traffic density, so that the speed law takes 
the functional form 

V{p) ^max (1 - p) (1 - p * 7/) . 
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The coefficient a in (|3.1I) is chosen so that j^rj — 1. The parameters a and h are the horizon 
of each driver, in the sense that a driver situated at x adjusts his speed according to the 
average vehicular density he sees on the interval [x — h^x — a]. To emphasize their roles, we 
select below the two situations 

a = -1/4 a = . . 

6 = 6 = 1/4. ^^-'^ 

In the former case, drivers look forward, while in the latter they look backward. We consider 
the initial datum 

13 3 

r \ ) 2 [-2.8, -1.8]^ ^ 4 ^[-1.2, -0.2]^ ^ 4^(0.6,1.0]^ ^ ^[1.5,+cx)[^ ^ ^ ^ 

representing three groups of vehicles lining up in a queue. 

The results in Section [2] ensure that for any G L-'^(R; [0, 1]), the Cauchy problem consist- 
ing of fll.ip - fl3.ip with initial datum admits a unique solution p — p{t^ x) attaining values 
in [0,1]. However, the qualitative behaviors of the solutions are rather different in the two 




Figure 1: Integration of (IIII])-(EI])-([331) in the two cases at times t = 0.05, 2.50, 5.01, 
7.50, 10.00. Above, drivers look backward while below they look forward. Note, already on 
the first column, the difference in the two evolutions, clearly due to the position of the support 
of 77. 

situations in (|3.2I) , see Figure [H Clearly, the evolution in the case of drivers looking forward 
(second line in Figured]) is far more reasonable, while the backward looking case leads to big 
oscillations in the vehicular density. 

3.2 Increase of the Total Variation and of the Norm 

This paragraph is devoted to show that Lemma 12.21 and the total variation bound fj2.8l) 
are, at least qualitatively, optimal. Moreover, the example below shows that the nonlocal 
equation (|l.ip does not enjoy two standard properties typical of ID scalar conservation laws, 
namely the maximum principle [1^ (iv) in Theorem 6.3], see also [12^ Theorem 3], and the 
diminishing of the total variation p4 Theorem 6.1]. 

In Remark 12.31 the assumption that f(p) = can not be replaced by v(p) = to ensure 
that the solution remains bounded between and p. Let p = 1 and consider fll.ip in the case 

IT /O 

f{p) = p, v{r) = l-r and r]{x) = a {{x - a){b - x)) X,.Xx) (3.4) 
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with a chosen so that j^r] — 1. Then, clearly, the initial data p{x) = 1 and p{x) = are 
stationary solutions to (jl.ll) - (j3.4l) . However, as the numerical integration below shows, the 
initial datum 

/(x) = 0.25y, .(^) + Xr ix) + 0.75Y, (x) (3.5) 

^ ^ ^ ^[-1.35,-0,95]^ ^ ^[-0.85,-0.25]^ ^ ^[-0.15, .25] ^ ^ ^ ^ 

which satisfies p^{x) G [0, 1] for all x G M, yields a solution p = p{t^x) that exceeds p = 1, 
showing that f|l.ip - fl3.4p does not satisfy the Maximum Principle, see Figure [21 




Figure 2: Numerical integration of (jl.ll) - fj3.4l) with a = 0, 6 = 0.2 and with the initial 
datum fl3.5l) at times t — 0.00, 0.25, 0.50, 0.75, 1.00. Note the sharp increase in both the 
L^-norm and in the total variation. 

Remark that the choice (|3.4I) leads to a flow in which is independent both from t 

and X. In the standard case of local scalar conservation laws, [1, Theorem 6.1] ensures that 
the total variation of the solution may not increase in time. The numerical integration below 
shows that the total variation of the solution to (|1.1I) - (|3.4I) may well sharply increase in a 
very short time, coherently with (|2.8I) . 

It is of interest to note that this behavior depends from the geometry of the support of 77. 
Indeed, a translation of the convolution kernel leads to very different solutions, see Figure [3l 
When the support is contained in R^, there is a sharp increase in the total variation. In the 




Figure 3: Total variation of the solution to (|l.ll) - (|3.4l) - (|3.5p versus time, in the three cases 
a = 0, 6 = 0.2; a = —0.1, h = 0.1 and a = —0.2, 6 = 0. Remark that the vertical scales in the 
leftmost diagram differs from that used in the middle and on the right. Indeed, the initial 
total variation is the same, 4, in all cases. 

other two cases, when sptTy is centered about the origin or contained in M~, there is a small 
increase in TV (p) for a small time interval, with a subsequent decrease. 

3.3 The Nonlocal to Local Limit 

In this section we use the algorithm (12. 5p to investigate the limit in which 77 tends to a Dirac 5, 
so that the nonlocal equation (II. ip tends, at least formally, to the local conservation law (II. 2p . 
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To our knowledge, no analytical result is at present available on this limit. 

Consider (|l.ip with the initial datum and the parameters below, see also Figure [H left: 



(3.6) 



f{t,x,p)^p v{r) ^ {1 - x^_^^^^{r) 

Va{x)^aa{a^-X^f/^X^_^^^^{X) Po{x) = !X[_i.8__i.3[(^)+X[_i.3^_o.8](^), 

and with the following choices for a: 

a = 0.25, 0.1, 0.05, 

with aa computed so that J^r]a{x)dx = 1. As limit case, we consider the standard conser 



Fundamental Diagram 



-2.0 -1.5 -1.0 -0.5 0.0 



0.5 1.0 1.5 




D.O 0.2 0.4 0.6 



0.8 1.0 



Figure 4: Left, the initial datum in (|3.6I) and, right, the flow (|3.6p used in (II. 2p . Note its 
change of convexity. 

vation law (II. 2p with / and v as in (|3.6I) , see also Figure [H right. In the integration below, 
the solution p attains positive values, so that after an easy modification of v on R~ we can 
assume that (12.21) holds. 

The resulting numerical integrations, carried out satisfying the CFL condition (12.41) , give 
the diagrams in Figure [5l In the limit case of (jl.2p , the chosen initial datum leads to the 
formation of a rarefaction wave, a shock and a mixed wave, due to the change of convexity 
of the flow, see the lowest line in Figure [5l The numerical integrations shown in Figure [5] 
qualitatively suggest that in the limit a ^ the solution to (Il.ip - (l3.6p converges to that 
of (I1.2p - (l3.6p . A more quantitative hint in this direction is in Figure [6l Using the algorithm 
above, we computed the solution pa to (Il.ip - (l3.6p for different values of a and the solution p 
to (I1.2p - (l3.6p . all at time t = 0.500. Figure [6] presents the plot of the L-'^-distance \\pa — p||li 
versus 1/a, see also table (13. 7p . 



1/a 

\Pa - pIIli 

1/a 

\pa - phi 

1/a 

\pa - pIIli 



4 5 6 7 

0.11166027 0.09569174 0.08373053 0.0743367 

9 10 20 40 

0.06049835 0.05526474 0.0282456 0.0122137 



0.06674645 
60 

0.0073903 



(3.7) 



80 100 150 200 250 

0.0054428 0.00449935 0.00354348 0.00319927 0.0030382 
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Figure 5: Integration of fll.ip - fl3.6p , first row with a = 0.25, second with a = 0.1, third with 
a = 0.05. On the last row, integration of fll.2p - fl3.6l) . The four columns display the times 
t = 0.5, 1.0, 1.5 and 2.0. The mixed waves are due to the non-convex flow (|3.6p . see Figure HI 



Convergence 




Figure 6: L-'^-distance between the solution pa to (|l.ip - (l3.6p for the values of a in (|3.7I) and 
the solution p to ([L2l)-(!3S1) at time t = 0.500. 
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4 Technical Details 

For any a, 6 G M, we denote I{a,b) = ]ci,b[ U ]b,a[. We use below the following classical 
notations: 

D^aj — ttj+i — ttj, — ~ -^^^j ~ ^j+i ~ + cij-i — {D^ — D~)aj 

and recall the trivial identities 

D^i^^jbj) = {D+aj)bj+i + {D+bj)aj, D-{ajhj) = {D-aj)hj + 

D'^{ajbj) = {D'^aj)bj + {D%j)aj + D+ajD+bj + D'ajD-bj. 

For later use, we note that the algorithm (|2.5p can then be rewritten as 



Proof of Lemma 12.21 Note that, by (|2.5I) . standard computations lead to 



p 



^+1 ^ (1 _ «n _ ^n)^n + ^n^n_^ + _ ^ (fj^./^lp^^ p]) " fj^l/2(P^ P,")) • (4-1) 



where 

^ ^^i^ — ^ 



Hi 



We now show that under condition (j2.4l) , the following inequalities hold: 



(4.2) 



«i e [0,1/3] 

/3| e [0,1/3] and A (f^+i/^l^.^ P") " fT-i/alPi , Pi )) < 3 ■ (4-3) 



G [1/3,1] 



Indeed, 



-\ ^p/(^'',^j+l/2,Cj+l/2) ^(c^'+i/s) + ^, 



So that 



"i < U^+a||viIlc 



> 



< 



(4.4) 
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Entirely similar computations lead to analogous estimates for /3". The bounds on 1 — a" — 
follow. The last term in ()4.3p . using (|2.6p and (|2.3p . is estimated as follows 



< 
< 
< h 



fit"", P]) v{c]^l/2) - fit""^ ^j-1/2, Pi ) 



+ 2||5p/| 



< C\\v\ 



< 



1 + 2 



r^o.h + 2\\dpfl 

II^p/IIl-) II^IIl- 







) 





- 3A 



Pj 



Using the bounds (14. 3p in (14. ip , we obtain 



proving the positivity of the discrete solution. 



□ 



Proof of Lemma 12.41 Thanks to the positivity of the discrete solution, it is sufficient to 
compute 



J 

J2hp]-h\(i\m ff+i/2(p^p^+i)-.iim fr_i/2(pr-i,pr; 



= IIp^IIli 

completing the proof. 

Proof of Lemma 12.51 For later use, estimate the quantity 

^ X] ^ \pk+l/2iVk-{j^l/2) - Vk-{j-l/2)) 



□ 



^■+1/2 S-1/2 



= ^IIP°llLi|h1lL^ 



Xk-j + 1/2 
Xk-j-1/2 

LOO 



77^(5)1 ds 



where Lemma [2^ was used. Using the same estimates as in the proof of Lemma [221 equal- 
ity (iJ]) yields 
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+x\f{r,Xj_y2,P])\\D-v{c]^y,) 

< (1 - a] - IIp^IIloo + a] \\p^^^ + \\p^^^ 



+A h 



7t _ r^f^ 



< IIp^IIloo + t 



provided 



(c I^IIloo + II^p/IIloo Ih'ILc IIpIIli Ih'Loo) Hp" 



IIp°IIli IIVIIl- • 



IL°° 



£ = C||HIloo + ||5p/||loo II^^'I 

A standard iterative argument completes the proof. 

Proof of Proposition 12.61 First, we write (12. 5p for j and for j + 1, subtract and get 

Pj+1 Pj - Pj+1 Pj 

Now add and subtract fjY3/2(/^j ' i^j+i) + ?)+i/2(^j-i' )' ^'^^^ rearrange to obtain 



(4.5) 
□ 



P]tl-P]^' = A]-XB] 



(4.6) 



where 



~^ (^1+3/2 (/^j+l ' /^j+2) ~ fj+l/2(/^j ' Pj+l) ~ fjV3/2(/^j ' /^j+l) + f)+l/2(/^j-l' /^j ) 
= f^3/2(P^Pi+l)-f^l/2(P^Pi^ (4.7) 

Consider first the term Aj. Recah fl4.2p and observe that, after suitable rearrangements, 
A] 

= I (Phi -Pi) 
+ {p]+2-Phi) 

+ (pI-pU) 

2 



' 1 A /(t",a;j+3/2,P")^(c ^+3/2) - Pj-l) ^(c"+i/2) ^ 

6^2 



^ - P]) 

+ {p]+2-PHi) 

+ ip]-p]-l 



6 ~ 2 



1 A fit"", Xj+3/2, p]) v{c]^s/2) - fi*""^ Xj+3/2,p]_i) ^(c^+3/2) 
6^2 - P"_i 



(4.8) 
(4.9) 
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+ 



(4.10) 
(4.11) 



Remark that the term (14. 8p equals defined in (14. 2p . Hence, it can be bounded 

using (14. 3p as follows: 



1 

< - 

- 6 



P]+2 - P]+l 



The estimate for the term (|4.9p is exactly as that of a" in Lemma [221 so that 



1 

< - 

- 6 



Pi - Pj-i 



Consider now the terms (|4.10l) and (14. lip : 

A / 

fPl+i ( \ 



dr 



2 



So that, passing to the absolute value 



1 

< - 
- 2 



LOO 



II- + II^p/IIl- II^IIl- I^IIlc 



n _ n 
Pj+1 Pj-1 



Grouping the estimates obtained we get 



An 



< 



1 + 



f 



WW-Loo \\p°\\li 



(4.12) 



B 



Pj+i-Pl 

We now turn to the term Bj' in (14. 7p . Since 



/(r , X^-+3/2, p^+i) v{Cj^3/2) - /(t^, X^- + l/2, P^+i) ^(Cj-+l/2) 

2 
2 



(4.13) 
(4.14) 
(4.15) 
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we 



consider the various terms separately. 

{j4J^ = V[Cj+l/2) ^ 

,f(,n ^ r,. - MCj+1/2) + v{Cj-l/2) 
+J{t ,Xj_^_i/2, Pj ) 



^ 2 



2 



where 



2 

2 2 



< 

- 2 



9x/r,o+i,P^)-9x/r,Ci-i,P,"; 



< 2 c/ftx,x,p,") 



dx 



< Ch^ 

where ()2.ip was used to get to the last hne. Moreover, 

2 

'^(cj+3/2) - vicj+1/2) _ vicj+1/2) + vicj-1/2) 
9 9 



= I {v'{Cj){c]^y2 - C^l/2) - ^'(Ci+l)(c^3/2 - C^l/2)) 

= ^ - ^''(O+l)) (c^+1/2 - cU/2) - ^^'(0+l)(c^3/2 + 2c^l/2 - C^l/2) 

= ^ ^"(e,) (0 - 0+1) (c^+1/2 - C^l/2) - ^ ^'(O+l) (c^3/2 + 2c^l/2 - C^l/2) • 

+ c^i 1 /o — 1 /o , and so using Young's 



Note that we have \(j — Cj+i\ ^ ~ 
inequahty, 



^(Cj+3/2) - 2^(Cj+i/2) + v{Cj-l/2) 

2 



2 1 

,n _ n \ _ — 



+ ^lh'llLc 



S+3/2 ^ ^S+1/2 S-1/2 • 

We now estimate the terms involving the discrete derivatives of in the expression above, 



14 



exploiting the regularity (12. 2p of 77. By (|2.7I) , we have 



S+1/2 ^--1/2 



5Z^Pfc+l/2 (^/c-(j+l/2) - ^A;-(j-l/2)) 



< 



< 



E^|^^J-l/2 



(4.16) 



Similarly, 



+ 2c^+l/2 ~ ^j-1/2 ^ /^fe 



-J-1/2 



i-1/2 



Vk-1/2 - 2?7fc+l/2 + '?fc+3/2 

v'{Ck+i)-v'{Ck)\ 



I /"U+I 



2/^'l|P"llLlh"||roo, 



(4.17) 



to complete the estimate of ()4.13|) we use the results above to bound the remaining terms 



/(t", P]) - fit"", P]) 



v{Cj+3/2) - v{cj+i/2) 



f{t'',Xj+i/2,p])-f{t^,Xj_y2,P]) 



< -hC 



< o^lk'llT.oollVllroollpl 



< -hC 



Li 



v{Cj+l/2) - v{Cj_i/2) 



< -h\\v'\ 
- 2 II I 



LOO IIVIIloo IIp^IIli 



and we are now able to complete the estimate of (|4.13p : 



(TO < h^C\\v\\^^ 



+ \\dpf 



pj IIloo 



Li 



c\\v\\^^ + {\\dpf\ 



LOO 



IL°^ W ' ML 



LOO 
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c \ 



Ili 



< \\v\\-^2,oo C + 



^pJ IIl- + ^ Ihllwi'- II/^'^IIli j 
We now pass to estimate fl4.14p and fl4.15p : 



gUD + 635]) = ^ (/(r, X,+3/2, P^+l) - , ^i+l/2, P^l)) ^(c,+3/2) 

Pi-l) - a;j--l/2, p]-l)) v{Cj+l/2) 

= \h (a,/(^^c,+l,P^+l)^(ci+3/2)-9x/(^^e„P^l)^'^ 

+^ (^/(^", a;j+i/2, ^''(7i+l) (ci+3/2 - Cj+1/2) 

-/(r, a:j_i/2, p]_i) v'{jj) {cj+i/2 - Cj-i/2 



for suitable G 



^j-l/2i Xj+1/2 



and 7j G i"(cj_i/2, Cj+1/2). Introducing ij G 

e Iip]_i,p]+i), 7j e /(7i,7i+i), 



Cj e I{p]-vP]+i)^ 7i e -f(7i,7i+i), 4' e J 3^7-1/2, 3:^+1/2 
5j e /(cj+3/2 - Cj+i/2,Cj+i/2 - Cj_i/2) and using (|2II]), (|4.16p, (|4.17p 

(I1I3D + dnsi) 



1 

+ 2 



^(7j+i)/i+ ^(7,+i) 



— 

+ 1 Pj-1 



+ 



Cj+3/2 - Cj+1/2 



K(7i)||5i||7i+i-7i| 



+ 



4' 0) \v'{jj)\ Cj+3/2 - 2Cj+i/2 - Cj_i/2 



< -hiCWv 



h + 



d"^ f 



Pj+1 Pj-1 



+ c\\v'\L \\po\\^, y 



+\ (clklLco IIpIIli IIVIIloo + II^p/IIloo II^'IIloo IIpIIli IIVIIlo 







LOO 



+2 



II Loo II J^oo T.o 



-C(l + ||r?'||j^^ IIpIIli 



+^ (c + II^p/IIloo II/IIli Ih'L- + 2 ll^p/ILoo) Ih'llwi.oo Ih'llwi.oo IIpIIli 







16 



1 

+ 2 



^2 f 

^pxJ 



h. 



I^IIl- p]+i p]-i 

The above bound allows to obtain the estimate for Bj: 

< lkllw^,oofc+(||9,/||j^^ + |)||ry'||^,^||p°||Li 



Pi 



+^c(l+||r7'||Loo IIpIIli) Ibllwi.oo |0, 



n 



IIpIIli c 



+ 



1 



pj+1 pj-1 



h. 



so 



that 



E^i < lbllw2.c«(^C'+(||9,/||j^^ + ^)||r?'||^,,^||p°||Lij||p°^ 

+^C'(l + ||r7'||Loo IIpIIli) Ibllwi.- IIpIIli^ 
+^ + IIV|Il- (2+IIp°IIli IIVIIloo)) Ikllwi- ll^'llwi,- 



LOO 

Recall now (jiSl) and (j3321) to obtain 



-Pi 



+ /C2 T 



where 

/C2 



1 II^p/IIl- II^'IIl- IIpIIli II^IIloo + 9px/ 



(4.18) 



2+l|p1lLih'|lL-))lh1 



l^'llwz.- IIp°IIli 



The estimate (|2.8p now foUows from standard iterative procedure. The proof of Proposi- 
tion ESI follows immediately. □ 



Proof of Lemma l2.7L We follow the same line as in [11, Section 3]. Using (|4.16p . Lemma 
Lemma 12.41 and Proposition 12.61 compute preliminarily 



1/2^ 
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\dj{t\ Ci, 0) ^(7i)| + , Ci, Ci) ^(7i) 

+ |/r,Ci,0)^'(7i) 
C||t;||Loop7/i+ ||9p/||loo I^IIl 

+ II^p/|Il°° Ih'll: 



Cj+1/2 S-1/2 

- 



I^IIl- max|pj,pj+i| 







Ili+2||9p/|U 


IIpIIli Ih'ILooll 


< C\\v\ 


II- IIpI 


Ili + 2||VILoo 





IL- iI^IIl- I ^2 ^ + X] - /^j 



The term Xlj^z 



admits an analogous estimate. Moreover, 



Using the above estimates and f|2.5l) we get 



,n+l _ n 



r 



D^[fir^^3-V2.Pl)v{d;_y,) 

< C{t)T 



jG2 



where 



at) = c\\v\\.^^\\p^^^ + 2\\d,f\L\\p"\\iy^ 



iLoo II ^ IlLoo 

PJ+1 - Pj 



completing the proof. 



(4.19) 



□ 



Proof of Proposition 12. 8L Fix n G N and for any sequence {p)jez define the transformation 
p !->• H{p) given by 



(4.20) 



H-{p) = p, - (f,'Vi/2(p„P,+i) - f]^i/2(P.-i, 
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where the functions fjYi/2 given by (|2.6I) , but where, instead of (|2.7I) , the sequence 
i^^^i/2)j^'^ is now an arbitrary fixed sequence. Thus, H^{p) depends only on Pj-i, Pj and 
Pj+i. Then, is monotone, in the sense that 



dH 



dp, 



^ >0, z = j-l,j,j + l. 



(4.21) 



The cases i = j =b 1 are easily verified. If i = j, using (|2.6p we find 
dm 1 A / 



dpj 



> 3 -^Fp/||l°o II^^IIl- 

> 



by the CFL condition ([231). The definition ()i:20|) of m and (f^ imply that for any k G 



p,- - fc| - A ( F,'Vi/2(Pi,Pi+i) - ^;-i/2(Pi-i,Pi) ) = ^7(P A fc) - if;(p V fc) 



(4.22) 



where k in the right-hand side above is understood as the sequence identically equal to k. 
The monotonicity condition (|4.21|) and the scheme ()2.5|) - (j2.6p ensure that 

H^{pAk)-H^{pVk) 

> H^ip)AH^ik)-H^{p)yH^{k) 

= sgn H^ip) -k + X (/(r, x,+i/2, k) - fit"", Xj_i/2, k) 

X H^{p) -k + \ (/(r, X,+i/2, k) - /(r, X,_i/2, k)) 

> sgn (if;(p) - fc) i/7(p)-fc + A(/(r,x,+i/2,fc)-/(r,x,_i/2,fc)) 
= i/7(p) - fc| + A sgn (if;(p) - k) [f{t\ x,-+i/2, fc) - f{t\ x,-_i/2, fc)) . (4.23) 

In the last inequality we used also the non-negativity of the function (a, b) i-)- (sgn(a -|- 6) — 
sgn(a))(a -|- 6). From and we conclude that 



(4.24) 



H^ip) -k -\pj-k\+\ (^FjYi/2(Pi'Pi+i) - Fj-i/2(Pj-i^Pj)) 

+Asgn(if-(p)-fe)(/(r,x,+i/2,fe)-/(r,x,_i/2,fe)) < 0. 



Consider now the numerical approximation p'J given by the algorithm (12. 5p . Then, we ap- 
ply (I4.24P to p^, with the sequence appearing in (I4.20p as given by the convolution (12. 7p . 
Observing that H^ip"") = we conclude that (I^TUll holds. □ 
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